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Abstract 

Accurate simulations of isotropic permanent magnets require to take the magnetiza¬ 
tion process into account and consider the anisotropic, nonlinear, and hysteretic material 
behaviour near the saturation configuration. An efficient method for the solution of the 
magnetostatic Maxwell equations including the description of isotropic permanent magnets 
is presented. The algorithm can easily be implemented on top of existing finite element 
methods and does not require a full characterization of the hysteresis of the magnetic ma¬ 
terial. Strayfield measurements of an isotropic permanent magnet and simulation results 
are in good agreement and highlight the importance of a proper description of the isotropic 
material. 
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1 Introduction 


Isotropic permanent magnets are used in a wide range of applications, like electric motors and 
actuators, speakers, as well as magnetic sensors [T|. Especially their lower production costs, 
and the additional capabilities due to the magnetizability in any direction, make them an 
interesting alternative to high-end anisotropic magnets. Due to the missing preferred direction 
the magnetization process has a strong influence on the final remanence magnetization within 
isotropic magnets. Accurate simulations of magnetic components therefore need to consider the 
non-linear, anisotropic, and hysteretic material behaviour at least close to the perfectly aligned 
saturation state. 

The most famous models for the macroscopic description of magnetic hysteresis are the 
Preisach model [2j as well as the Jiles-Atherton model f3]j. Extensions and combinations of 
those models have been proposed, where the applicability to the vector hysteresis for isotropic 
magnets has been demonstrated @10 0 E]. The integration of those material descriptions 
into finite element simulations is based on either Fixed-Point iteration or the Newton-Raphson 
method for the solution of the resulting nonlinear systems of equations umm- 

Although those methods are well suited for the simulation of dynamic hysteresis effects, the 
complete description of the hysteresis involves a large computational effort and is not always 
needed. Considering permanent magnetic materials, the correct description of the hysteretic 
behaviour near the saturation configuration is often sufficient. In addition to the higher com¬ 
putational effort, the complete characterization of a magnetic material requires complicated 
reversal curve measurements, which are not always available for commercial magnets. 

Within this work a new method for the description of hysteretic behaviour of isotropic 
permanent magnets is introduced. The presented algorithm allows to describe small deviations 
from the saturation configuration in a phenomenological way by means of a single parameter - 
the freezing field. The method is very efficient and easy to implement on top of existing magnetic 
finite element simulations. Moreover there is no need for complete reversal curve measurements 
of the magnetic material. 

The structure of the paper is as follows: Section [2] summarizes the numerical method used 
to simulate the behaviour of anisotropic magnetic components. Within Section [3] a simple phe¬ 
nomenological algorithm is introduced which allows to extend the method to isotropic magnets. 
Finally simulations of the magnetic strayfield of an isotropic permanent magnet are compared 
with measurement data within Section HI 

2 Anisotropic Magnets 

Anisotropic magnets are produced with a magnetic field applied during the production process. 
The magnetic grains orient into the direction of the applied field. After the sintering or bonding 
process the applied field is turned off, but the prescribed orientation of the magnetic grains 
remains. Thus, the crystallites are fixed and even very high applied fields do not affect their 
orientation. There exists a global preferred direction and it is possible to define a homogeneous 
intrinsic anisotropic material law for the whole magnet. Assuming a uniaxial anisotropy requires 
the distinction of longitudinal and transversal behaviour. 

For sake of simplicity the global material law will be linearized around the remanence vector 
B r . Assuming the longitudinal easy axis is oriented into the 2 axis, results in the following 
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constitutive law 


B(H) = B r + /iH 

/ 0 \ (fi_L 0 0 

B r = I 0 1 /I = I 0 //j_ 0 

V-Br/ \ 0 0 / 1 || 

, with the magnetic flux density B, the total field H, the remanence vector B r and the perme¬ 
ability tensor /x. Nonlinear behaviour can be described by updating the linearized /i depending 
on the actual local field H. For typical applications the demagnetization field shifts the local 
operating point in a region of the B(H)-curve, which can be well described with a linear model. 

It is important to distinguish between the B r -vector, the magnetic flux-density without any 
total magnetic field, as well as the magnetic remanence B(H ext = 0), which requires a vanishing 
external field but includes the magnetic strayfield. Measurements typically yield the magnetic 
remanence, unless a closed loop measurement is performed, which eliminates the magnetic 
strayfield and therefor allows direct measurement of the B r -vector. The terms ”B r -vector” as 
well as ” magnetic remanence” are used in the following to distinguish these two quantities. 

The actual magnetization state of the permanent magnet can be calculated by means of 
the homogeneous, magnetostatic Maxwell equations. These are solved efficiently by FEMME 
Maxwell, which uses a modified version of the method presented in m- The (maybe inhomoge¬ 
neous) magnetic material is contained within the magnetic domain 12 and is treated by a finite 
element (FEM) approach. The homogeneous outer region is handled by a boundary element 
(BEM) approach which is able to deal with open-boundary problems very efficiently. Maxwell’s 
equations as well as the jump conditions at the boundary of the magnetic parts are given by 


rot H = j in [R 3 (2a) 

div B = 0 in R 3 (2b) 

[B] • n = 0 on d!2 (2c) 

[H] x n = 0 on d!2 (2d) 



, where [.] denotes the jump of a quantity over the boundary 912 of the magnetic domain, 
and n is the outward pointing surface unit normal. Introducing a reduced scalar potential u 
and expressing the magnetic field created by the current density j as H ext finally leads to the 
following system of equations 


H = Hext - Vu 


(3) 

V • (fi Vm) = V • (/j, Hext + B r ) 

in 12 

(4a) 

Au = 0 

in [R 3 \ 12 

(4b) 

N =0 

on d!2 

(4c) 

[Vm • fi ■ n] = (Hext • M + B r ) • n 

on 912 

(4d) 


3 Isotropic Magnets 

In contrast to anisotropic magnets described in the previous chapter, isotropic magnets are 
produced without an applied field. The magnetic grains are oriented randomly inside of the 
magnet, and thus the material has no preferred direction. An advantage of this type of magnet 
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is that it is easier to produce and therefore less expensive. Furthermore those magnets can be 
magnetized into any desired direction due to their isotropic behaviour. On the other hand the 
magnetic remanence, the coercive field as well as the energy product of isotropic magnets are 
significantly reduced compared to anisotropic magnets of the same material. 

One additional drawback of isotropic magnets, which should be emphasized in this chapter, is 
that there exists no intrinsic anisotropy axis of the material. Instead the anisotropic behaviour of 
the magnet is created by the application of a magnetic field pulse after the production process. 
For large enough fields the average magnetic moment within the magnet will point into the 
direction of the spatially varying total field. As a consequence the local B r -vector within the 
magnet is not only modified at the peak of the magnetization pulse, where the total field is 
almost equal to the applied field, but also during the decay of the applied field pulse. The 
demagnetization field gets more and more influence and the B r -vector state gets increasingly 
inhomogeneous during the reduction of the applied field. When the total field at a particular 
location comes below a certain threshold, the B r -vector freezes. 

This effect leads to a shape dependent B r -vector state and therefore requires taking into 
account the specific magnetization process. This can be realized by introducing a critical freezing 
field -ff freeze and modify the B r -vector locally only if the freezing condition |H| > i7f reeze is 
fulfilled. 

-of | B r | H if |H| > freeze ( r\ 

r \ B r otherwise 

Assuming a simple Stoner-Wohlfarth model for the individual magnetic grains indicates a 
freezing field between H c /2 < E/f reeze < H c , where H c is the is the coercive field of the individual 
particle. Starting from a perfectly aligned B r -vector state the applied field is gradually reduced 
and the total field is calculated from equations © and 0- The field steps need to be chosen 
small enough to create self-consistent results. In the following chapter simulation results using 
the presented procedure are compared with measurement data. 

4 Comparison with Measurements 

The importance of proper simulation of the magnetization process should be demonstrated 
by a strayfield measurement of an isotropic magnet as it is used as bias magnet in magnetic 
sensor applications. The dimensions of the magnet used for the measurement is approximately 
7.0mm x4.5mm x5.5mm. The commercially available plastic bonded magnet consists of a 
NdFeB-composit material specified in Tbl. [l] and is directly moulded onto an integrated sensor 

m- 


Rare earth magnets 
(injection moulded) 

Energy 
product 
(B * H) max 
[kJ/m 3 ] 

Rem 

b y 

[mT] 

anence 

Hf 

[kA/m] 

Coercivity 

H C M 

[kA/m] 

Relative 
permeability 
M||/Mo 
[1] 

Curie 

Temperature 

T c 

[°C] 

isotropic NdFeB 

43.0 

505 

340 

800 

1.15 

305 


Table 1: Summary of the magnetic properties of the plastic bonded rare earth magnet material used 
for the comparison of measurements and simulations. 

For the measurement the magnet is magnetized using an impulse magnetizer, which creates a 
peak magnetic field of /io7f max = 4T. After magnetization in z direction the magnetic strayfield 
along the rr-axis is measured 1.0 mm above the magnet (see Fig. [5]) using a Hall probe. 

Three different simulations are performed and compared with the measurement results: 
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(a) Anisotropic magnet model using FEMME Maxwell (11] 

A homogeneous remanence vector B r = 0.51 T is assumed within the whole magnetic region. 
A homogeneous, isotropic fi = = /iy = 1.12 /jLq is chosen since only the longitudinal 

fi\\ is known from the magnet datasheet. Additional simulations using an anisotropic /i, 
with ji l = 1.19 fi§ deduced from measurements, show that the effect of the anisotropic 
permeability is negligible compared to the shown remanence-freezing effect, which justifies 
the choice of an isotropic permeability. The magnetization state of the magnet is calculated 
without external field by means of equations 0 combined with the material law 0 - 

(b) Anisotropic magnet model using ANSYS Maxwell 3D (13j 

The commercial and well known electromagnetic field simulation software ANSYS is used to 
validate the results of the previous simulation. It allows to define anisotropic, hard-magnetic 
materials using the same parameters M r and i± as above. 

(c) Isotropic magnet model using FEMME Maxwell with remanence-freezing 

For the simulation of the isotropic model it is necessary to incorporate the magnetization 
process. Therefore an external field of /igi^ext = Mo^max = 4T is applied and successively 
reduced in 0.1 T steps until zero external field is reached. Within each step an anisotropic 
model with the same parameters as above is solve. After each field-step the freezing condi¬ 
tion 0 is checked and the local B r -vector is updated for each element. 

A comparison of the anisotropic and the isotropic simulations is shown in Fig. [l] and [2j 
respectively. In contrast to the homogeneous B r -vector assumed for the anisotropic simulation 
the isotropic B r -vector state gradually degrades into a flower state while the applied field 
is reduced. When the external field is reduced below i7 ex t ~ ^freeze the B r -vector state is 
completely frozen within the whole magnet. Choosing Hf reeze = 0.6 T yields good agreements 
with measurement data for various geometries. Due to the shape of the magnetic strayfield the 
remanence-freezing effect becomes visible mostly on sharp edges of the geometry. 

Finally the comparison of the magnetic strayfield above the magnetized magnet is shown 
in Fig. [3j The anisotropic models simulated with FEMME Q and ANSYS coincide very 
well, which validates the implementation of the anisotropic model presented in Chapter [2] within 
FEMME. Furthermore it can be seen that the measurement data shows significant deviations 
from the anisotropic behaviour. The simulations using the isotropic magnet model ([c]) presented 
in Chapter [3j leads to a much better description of the measurement results. The derived flower 
state for the B r -vector state amplifies the influence of the magnetic strayfield onto the actual 
remanent magnetization. 

5 Conclusion 

An efficient method for the description of isotropic magnetic materials is presented. Comparison 
with measurement data reveals that in contrast to anisotropic magnets the magnetization pro¬ 
cess has a significant influence on the final remanent magnetization of isotropic magnets. This 
effect can approximately be described by the introduced remanence-freezing method, which 
leads to a very good agreement with the measurement results. The method can easily be 
integrated into existing anisotropic magnetic simulations. Additionally no detailed hysteresis 
measurements of the magnetic material are required, which simplifies the utilization for practical 
applications. 
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Figure 1: Simulation results using an anisotropic magnet model ([a]). A homogeneous B r -vector is 
assumed within the whole magnetic region (top). The resulting magnetic polarization (bottom) shows a 
weak flower state, but all magnetic moments are mostly aligned with the prescribed preferred direction. 
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